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We perform an in-depth study for mean first-passage time (MFPT) — a primary quantity for random walks 
with numerous applications — of maximal-entropy random walks (MERW) performed in complex 
networks. For MERW in a general network, we derive an explicit expression of MFPT in terms of the 
eigenvalues and eigenvectors of the adjacency matrix associated with the network. For MERW in 
uncorrelated networks, we also provide a theoretical formula of MFPT at the mean-field level, based on 
which we further evaluate the dominant scalings of MFPT to different targets for MERW in uncorrelated 
scale-free networks, and compare the results with those corresponding to traditional unbiased random 
walks (TURW). We show that the MFPT to a hub node is much lower for MERW than for TURW. However, 
when the destination is a node with the least degree or a uniformly chosen node, the MFPT is higher for 
MERW than for TURW. Since MFPT to a uniformly chosen node measures real efficiency of search in 
networks, our work provides insight into general searching process in complex networks. 

Random walks in complex networks have been heavily studied in the past years' ' due to their wide range of 
applications in science and engineering''. Recently, continuously increasing efforts have been devoted to 
maximal-entropy random walks (MERW)^ also called Ruelle-Bowens random walks'"'", where all walk- 
ing trajectories from given starting and ending points of a given length are equiprobable. In this sense, MERW is 
the most random of random walks, which maximizes the entropy rate'^ ''' and is in striking contrast with the 
traditional unbiased random walks (TURW) and other biased random walks. The unique diffusion process of 
MERW leads to several unusual phenomena, such as localization of stationary distribution"^ and fast relaxation'''. 

As a powerful tool, MERW has been fruitfully applied in various fields. For instance, the localization phe- 
nomenon of stationary distribution for MERW makes it a good measure of centrality that is more discriminating 
than some classic centrality measures, e.g. PageRank, in the sense that it can distinguish evidently those nodes that 
PageRank regards as almost equally important'^"'. Furthermore, since MERW incorporates both network struc- 
ture and eigenvector centrality of nodes, it was also applied to establish a new algorithm of link prediction, which 
outperforms various supervised and unsupervised techniques of link prediction, on most test databases'^. In 
addition, MERW has also found applications in optimal sampUng algorithm'", demographic stability of popu- 
lation", community detection""". 

A fundamental quantity related to random walks is first-passage time (FPT)^'"^'', which is the expected time 
required for a random walker starting from a source point to a given target point. The mean first-passage time 
(MFPT) is defined as the average of FPTs over all source nodes in the network, which is a useful tool to analyze the 
behavior of random walks. The importance of MFPT originates from the essential role played by first encounter 
features appearing in various real situations, such as lighting harvesting^^"^' and target search^"'^^. The MFPT can 
also serve as a significant indicator measuring node importance™ and efficiency of trapping process"". It is thus of 
utmost importance to study MFPT for different random-walk processes. Thus far, MFPT has been intensively 
studied for TURW^^ '"' and some biased random walks'" ''*, while related research about MFPT for MERW is stiU 
much less, although the particular diffusion process of MERW is suggested to significantly affect the leading 
behavior of MFPT. 

In this paper, we propose a theoretical framework for MERW in complex networks and perform an in-depth 
study on the MFPT for MERW to a given target. We derive an explicit expression of FPT for MERW from one 
node to another in any connected network in terms of the eigenvalues and eigenvectors of adjacency matrix for 
the network. Based on the obtained representation for FPT, we further deduce an exact formula for MFPT to an 
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arbitrary target node. Moreover, for uncorrelated networks, we also 
provide an analytical expression of MFPT for MERW at the mean- 
field level, using which we obtain the leading scalings of MFPT for 
uncorrelated scale-free networks with various degree exponent y, and 
show how the MFPT scales with the network size. 

For MERW in uncorrelated scale-free networks, we study the 
MFPT for three representative cases with the target node being 
located at a hub node, a node with the smallest degree, and a node 
uniformly chosen from the system, respectively. For all the three 
cases, we derive analytically the leading scalings for MFPT, all of 
which depend on the degree exponent y that characterizes the het- 
erogeneity extent of scale-free networks. Our results indicate that for 
the last two cases that the target is placed at a smallest degree node or 
a uniformly selected node, the leading scalings resemble each other, 
but both scalings are considerably larger than that corresponding to 
the first case when a hub is the target. 

We also compare the obtained results of MFPT for uncorrelated 
scale-free networks with those corresponding to TURW. We show 
that when the target is frxed at a hub node, the MFPT for MERW is 
much less than that for TURW. On the contrary, when the target is 
placed at a smallest degree node or a randomly chosen node, the 
MFPT for MERW is larger than that associated with TURW. 
Therefore, in comparison with TURW, the special diffusion process 
of MERW has a stronger influence on the efficiency for searching a 
target in heterogeneous networks, making the process considerably 
more efficient for finding hub node but less efficient for locating a 
node with small degree or a randomly chosen node. 

Results 

Explicit expressions of MFPT for MERW. Throughout the paper, 
the random walk processes considered are defined in a connected 
undirected graph G with N nodes and E edges. The connectivity of 
nodes is described by the adjacency matrfac A, in which the entry = 
1 if nodes ; andj are adjacent, and Ujj = 0 otherwise. Then, the degree 

of a node ! is fcj = 2__^^i Si)- Let /'.i,22, ■ ■ ■ j/'-n be the N eigenvalues of 

A, rearranged as li> l2> ■ ■ ■ > ^n, and let /(j, fi2, ■ ■ ■ , fi/^he their 
corresponding mutually orthogonal eigenvectors of unit length, 
where /j, = (/.(,i,/.(,2, • • • Then, matrix A can be decomposed 

as the following spectral form: 

A = Udiag[Ai,;.2,---,AN]UT, (1) 

where V = {fii,iJ.2, ■ ■ ■ ,A'n) is an orthogonal matrix, obeying 

UU^ = U^U = I, (2) 

where I is the identity matrix. 

Using above notations, we can introduce the MERW that is char- 
acterized by a unique stochastic matrix P, the ijth entry of which is 
given by 



^-1 A'l, ' 



where Ai is the principal eigenvalue of matrix A, and f.iij is the ith 
entry of the principal eigenvector corresponding to Aj. This guar- 
antees that MERW maximizes the entropy of a set of trajectories with 
a given length and end-points and leads to the maximal entropy rate 
of such processes\ The stationary distribution of MERW is" 



71= (7ti,7i:2, ■ ■ ■ TInV = (/^ll'/'l2> ■ 



(4) 



The MERW is biased, which is different from TURW, where the 
transition probability = fly//c, from a node i to one of its neighbor- 
ing nodes j is identical. 

The main quantity we are interested in the paper is MFPT. Notice 
that MERW in an arbitrary connected binary network can be repre- 
sented as generic random walk in a corresponding weighted net- 



work*^. The !;'th element of the generalized adjacency matrix 
(weight matrix) W for the weighted network is defined by Wy = 
ajji^iijUy, which specifies the weight of the edge connecting nodes 
and j. In this weighted network, the strength""" of a node i is given by 

Sj = ^^.^ = and the total strength of all nodes is 

s= 2—i - \ ^' "''■i' ^'^^ generic random walks in this weighted net- 
work, the transition probability is defined as 



^1/^1,- 



(5) 



which is equal to transition probability, given by equation (3), for 
MERW in the original graph. The equivalence between the two ran- 
dom walks allows to determine the MFPT for MERW in a graph by 
evaluating the corresponding quantity for generic random walks in a 
related weighted network. 

For MERW in a network, let Ty denote the FPT from node ; to 
node j. Without loss of generality, let; be the target node, and let Tj be 
the MFPT to node j. Then, by definition, the MFPT Tj is given by 



' N 



1 



(6) 



Based on the equivalence between MERW and corresponding gen- 
eric random walks, we derive an exact expression (see Methods) for 
Tjj in terms of the eigenvalues and eigenvectors for adjacency matrix 
A: 



T 
^ 11 



1 Ai 



(7) 



Plugging equation (7) into equation (6), we arrive at an expression of 
MFPT Tj for MERW in a general graph with the deep trap fixed at an 
arbitrary node j, given by 



l-'-ki 



(8) 



Equation (8) provides a universal formula of MFPT to any node for 
MERW in an arbitrary network. Although it involves computing 
eigenvalues and eigenvectors of adjacency matrix, which puts heavy 
demands on time and computation resources for large networks, it 
can be utilized to check the results for MFPT obtained by other 
approaches, at least for small networks. Besides, equation (8) can 
also be used to compute the exact average (T) of Tj over all JV targets: 



1 N ^ N 



N{N- 



(9) 



[hi 
Uhi, 



which is exactly the MFPT when the target is uniformly distributed. 

A drawback for equation (8) is that by using this spectral technique 
it seems very difficult, even impossible, to obtain the leading behavior 
of MFPT Tj characterizing the random- walk dynamic process. Thus, 
it is important to seek alternative techniques of evaluating MFPT Tj 
even for particular networks, which are not computationally 
demanding but are valid to estimate the scaling of MFPT. 
Fortunately, for uncorrelated networks, we can derive an expression 
of MFPT for MERW at the mean-field level, as well as its dominant 
scaling for scale-free networks. The details will be given below. 

Theoretical prediction of MFPT for MERW in uncorrelated 
networks. We now consider the MFPT for MERW in uncorrelated 
networks, where the degree-degree correlations between adjacent 
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nodes are absent. In a recent work*', we have shown that, for generic 
random walks in uncorrelated weighted networks, the MFPT to node 
j can be represented in terms of the strengths of the nodes as 



s 



(10) 



Plugging Sj = Aifiy and 5 - 
node; for MERW: 



Ai into equation (10) gives the MFPT to 



1 



(11) 



Thus, in order to obtain Tj, it is sufficient to determine ^ij. Although 
the evaluation of eigenvectors of a general matrix is very hard, for 
uncorrelated networks we can approximate f.iy at the mean-field 
level (see Methods) as follows: 



(12) 



Substituting equation (12) into equation (11), we reach a theoretical 
approximation of MFPT for MERW to node j: 



T,- 



kj 



(13) 



In Fig. 1, we report both the exact numerical results and theoretical 
approximate results of MFPT for MERW taking place in Erdos- 
Renyi (ER) network''" and Barabasi- Albert (BA) network*', with 
both results being generated by equations (8) and (13), 
respectively. Figure 1 shows that the theoretical predictions agree 
well with the numerical results. From Fig. 1, we can also find that 
for nodes sharing identical degree, MFPT for MERW distributes in a 
broader range in BA network than in ER network. This difference lies 
in the structure of the networks. Since BA networks is heterogeneous, 
the component of leading eigenvector localizes at hub nodes and 
their neighbors^. For nodes having the same degree but different 
neighbors, their MFPT differ widely. For example, for two leaf 
nodes in treelike BA network that are linked to a hub node and a 
small-degree node farther from the hub, respectively, the MFPT to 
the leaf connected to a hub is much less than the MFPT to the other 
leaf While for ER network, it is almost homogeneous, so the 
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Figure 1 | MFPT to a given node for MERW in ER network (a) and BA 
network (b). Black circles represent the numerical results obtained by 
equation (8), and each red triangle stands for the average of numerical 
values for Tj over all nodes having the same degree kj. Straight lines are the 
theoretical approximation generated according to equation (13). 



disparity for MFPT to different nodes with identical degree is 
relatively indiscernible. 

In order to better understand the behavior of MFPT in inhomo- 
geneous networks, in the sequel, grounded on the theoretical 
approximation in equation (13), we will analytically evaluate the 
leading scaling of MFPT for MERW in uncorrelated scale-free net- 
works, aiming to unveil the effects of target location on the MFPT for 
MERW, as well as the difference between MERW and TURW in 
terms of the MFPT. 

Leading scalings of MFPT for MERW in uncorrelated scale-free 
networks. Extensive empirical studies™ have shown that most real- 
world networks exhibit the striking scale-free property*', 
characterized by a power-law degree distribution P{k) ~ k^" with 
y > 2. In this section, we will study the leading scalings of MFPT for 
MERW in uncorrelated scale-free networks. We wiU examine the 
dominant scalings of MFPT for three representative target 
problems, with the target being a hub with the highest degree, a 
node with the lowest degree, or a node selected uniformly. Our 
goals are twofold. One is to uncover of the influence of target 
location or degree on the behavior of MFPT. The other is to find 
the scaling difference of MFPT between MERW and TURW. 

Scaling of MFPT to a huh node. Let fcmax denote the degree of a hub 
node, and Th the MFPT to this hub node. Then, by equation (13), 



/j2 



(14) 



The numerator in equation (14) can be evaluated as 

NP{k)k^dk 



( Nk^-' 



2<y<3, 
Nln/cmax, y = 3, 



(15) 



N, 7>3. 



Note that in a scale-free network with N nodes and power-law degree 
distribution P{k) ~ k^'', the largest connectivity can be esti- 
mated as^' 



(16) 



Combining equations (14-16), we can obtain the leading scaling of 



' N°, 2<y<3, 
Th^I InN, y = 3, 

, y>3. 



(17) 



Thus, the extent of inhomogeneity, characterized by the degree expo- 
nent y, of scale-free networks strongly affects on the MFPT Th to a 
hub node for MERW. For 2 < y < 3, Th is approximately equal to a 
constant; for y = 3, Th grows logarithmically with the network size N; 
while for y > 3, Th grows sublinearly with N. 

We have checked our approximate results for against numer- 
ical values obtained according to equation (8) for uncorrelated scale- 
free networks with various values of y, namely, y = 2.5, )' = 3 and y = 
3.5. The considered network with y = 3 is the BA model; while the 
networks with y = 2.5 and y = 3.5 are generalizations of the BA 
modeP^. The comparison for theoretical and numerical results is 
shown in Fig. 2, which indicates that for different values of y and 
network size N, the analytical predictions from equation (17) agree 
with those numerical results given by equation (8). It should be 
mentioned that for y = 2.5, the prediction is only valid for large 
networks, since in the process generating scale-free networks with 
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2 < y < 3, multiple and self-connections are forbidden, which could 
introduce degree correlations in resultant networks^''^"*, leading to a 
deviation of numerical values from theoretical prediction. 

We next show that the behavior of MFPT for MERW in scale-free 
networks is quite different from that for TURW in the same net- 
works. Similar to MFPT in weighted networks, the MFPT to a hub 
node for TURW in uncorrelated scale-free networks can be esti- 
mated by 



Th 



where the numerator kj can be approximated as 

NP{k)kdk~N. 



(18) 



(19) 



Considering equation (16), the leading scaling of for TURW is 

rH~N('-2)/'''-i', (20) 

which scales sublinearly with the network size Nbut decreases with )', 
a result consistent with that previously obtained^"" using a different 
approach. In Fig. 2, we also plot the approximation for in equa- 
tion (20) against their corresponding numerical values for TURW 
generated by the method in'^, both of which agree well with each 
other. 

Equations (17) and (20) show that the MFPT to a hub for MERW 
and TURW in uncorrelated scale-free networks display rich but 
distinct behavior. For both MERW and TURW, the MFPT depends 
on the exponent y: lower y corresponds to smaller MFPT. Moreover, 
in the whole range of y, the MFPT for MERW is smaller than its 
counterpart for TURW. The root of the difference of MFPT between 
TURW and MERW is attributed to their local transition probabil- 
ities. For TURW, the transition probability pij from a node i to a 
neighboring node j is identical, while for MERW, the transition 
probability ispy = fiy/iAiiiii) ~ kj/{).iki), proportional to the degree 
of the neighboring node j. Thus, a walker visits a hub node more 
quickly in MERW than in TURW. 

Scaling of MFPT to a node with the smallest degree. We now study the 
MFPT in uncorrelated scale-free networks when the target is a node 
with the smallest degree. According to equation (13), the MFPT to a 
smallest node can be represented as 



(21) 



where fcmin denotes the degree of a node with the least degree. For 
sparse scale-free networks, their average node degree is a small con- 
stant™. Hence the minimal degree k-^^ can be regarded as a smaller 
constant. Then, recalling equations (15) and (16), the dominant 
scaling of Ts can be approximated by 

N2/(;-i), 

(22) 



NlnN, 



2<j'<3, 

y=3, 



'>3, 

which is supported by extensive numerical simulations, see Fig. 3. 

Equation (22) indicates that the MFPT Tg for MERW in uncorre- 
lated scale-free networks also exhibits rich behavior relying on the 
degree exponent y. When 2 < )' < 3, Ts varies superlinearly with the 
network size N; when y = 3, Tg scales with Nas Nln N; when )' > 3, Tg 
behaves linearly with N. 

Although for both cases that the target is either at a hub node or at 
a smallest node, the MFPT is influenced by the degree parameter y, 
the dependence relation of MFPT on y is quite distinct, as can be seen 
from equations (17) and (22). Furthermore, in the whole range of 
2 < y < 00, Th is much smaller than Ts- As a result, for MERW in 
uncorrelated scale-free networks, generating all paths with identical 
probability is disadvantageous for the walker to explore nodes with 
small degree. 

In addition, the MFPT Tg for MERW is also different from that of 
TURW, where the MFPT Tg to a node with the smallest degree can be 
estimated as 



Y^N ■ N 



■N, 



(23) 



where equation (19) is used. This analytical solution is conformed by 
numerical results shown Fig. 3. 

Equations (23) implies that for TURW, the MFPT Ts scales line- 
arly with network size N, independent of y, which is totally different 
from the behavior of MFPT, provided by equation (22), correspond- 
ing to MERW. From equations (22) and (23), we can observe that for 
2 < y < 3, the MFPT Tg for MERW is larger than that for TURW; 
and that for y > 3, although the leading scaling of MFPT Tg grows 
linearly with network N for both MERW and TURW, the values of Ts 




100 



N 



1000 



Figure 2 | MFPT to a hub node as a function of the network size N. The 
filled symbols are the data of numerical results. The lines correspond to 
theoretical predictions provided by equations (17) or (20). 




Figure 3 | MFPT to the node with the least links. The filled symbols stand 
for the numerical data, each being an average of MFPT over all nodes 
having the smallest degree; while the lines refer to the theoretical 
approximations provided by equations (22) or (23). 
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for MERW is greater than those corresponding to TURW, which can 
be seen from Fig. 3. Thus, for the case of target node located at a node 
with the smallest degree, performing MERW is substantially slowly 
to arrive at the destination than performing TURW, which is in stark 
contrast with the the case when the target node is a hub node, for 
which, performing MERW is more efficient than TURW to find the 
target. This phenomenon can also be accounted for by local trans- 
ition probability from a node to one of its neighboring nodes with 
small degree, which is relatively smaller for MERW than for TURW. 

Scaling of MFPT to a uniformly chosen node. The above studied MFPT 
to a particular target node in a network is often not looked upon as a 
general dynamical property of the network""" '^. Instead, the average of 
MFPT over all targets reflects some global characteristics such as the 
efficiency of searching process. Thus, it is of significance to compute the 
target averaged MFPT. Next, we address random walks in an uncorre- 
lated scale-free network with the target being a uniformly selected node 
in the network. In such situation, the MFPT (T) is defined as the 
average of FPTs over all pairs of nodes in the network, which involves 
a double average: The former is over all the source nodes to a given 
target node, the latter is the average of the first one. That is. 



;=1 



(24) 



Next, we determine (T) for MERW and TURW, respectively. 
For MERW, plugging equation (13) into equation (24) gives 

, N , N 

;=1 J 1=1 

where the term kj'^ can be estimated by 



(25) 



N 



j=i ^ 



NP{k)k-^dk~N. 



(26) 



Considering equations (15), (16) and (26), the quantity (T) follows 

f]V2/()-i), 2<y<3, 
NlnN, y = 3, (27) 
N, y>3, 

which is consistent with the numerical results, see Fig. 4. 

By comparing equation (22) and (27), we observe that for MERW 
the MFPT (r) exhibits similar behavior as that of Ts. This phenom- 
enon can be heuristically understood from the structure of scale-free 
networks, where the fraction of nodes with small degrees is very high. 
Moreover, the MFPT to a small-degree node is much larger than that 
of large- degree node. Thus, (T) and Ts resemble in the leading beha- 
vior, which means that the dominant scaling of Ts to a small-degree 
node is representative of MERW in scale-free networks. 

We proceed to uncover the difference for (T) between MERW and 
TURW. For TURW, the quantity (T) can be approximated by 



where the term 



N 



j=i J 1=1 



can be estimated as 



NP{k)k-^dk^N. 



(28) 



Recalling equation (19), we have 

(r)~N, 



(29) 



(30) 
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Figure 4 | MFPT to a target node uniformly selected from the whole 
network. The filled symbols are the numerical results generated by 
equation (9); while the lines correspond to the theoretical predictions 
given by equations (27) or (30). 

a scaling similar to that of Ts for TURW. In Fig. 4, we plot the 
numerical results of (T) versus theoretical prediction in equation 
(30) for TURW in scale-free networks with different y, which are 
consistent with each other. 

On the other hand, for uncorrected networks, the relation 

l/Tj^l holds. Then, according to inequality of arithmetic 

and geometric means, we can deduce a lower bound of (T) for uncor- 
related networks: 




N 



(31) 



which provides a minimal scaling for (T) in uncorrelated networks. 

Equations (27) and (30) show that although the scaling of (T) for 
MERW and TURW behaves differently, the optimal linear scaling for 
(T) can be achieved both for TURW in scale-free networks with 
arbitrary 2 < y < =o and for MERW in scale-free networks with y 
> 3. In the end, we stress that although in the case of y > 3, for both 
MERW and TURW, (T) can reach the minimal scaling, the cofactor 
of the dominating scaling for (T) is larger for MERW than for 
TURW, which can be seen in Fig. 4. Therefore, in the whole range 
of 2 < y < =o, the value of <T) is higher for MERW than that for 
TURW. 

Discussion 

In summary, we have presented a comprehensive and systematical 
analysis of MFPT for MERW in complex networks. We have pro- 
vided an explicit expression of MFPT for MERW in a general net- 
work with a target node being located at an arbitrary node, which is 
provided in terms of eigenvalues and eigenvectors of the adjacency 
matrix for the network. Moreover, for MERW in an uncorrelated 
network, we have given an alternative theoretical prediction for 
MFPT at the mean-field level, which is devoid of calculating the 
eigenvalues and eigenvectors but gives good approximation for 
MFPT that are confirmed by extensive numerical results. 

Applying the mean-field approximation formula, we have further 
addressed the leading behavior of MFPT for MERW in uncorrelated 
scale-free networks with a given target and various degree exponent 
y, focusing on three representative cases with the target being a hub 
node, or a node with the least links, or a node chosen uniformly. For 
all the three cases, the MFPT is dependent on the degree of the target, 
as well as the degree exponent y. We have also performed a compar- 
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ison of the obtained results for MERW with those corresponding to 
TURW. For the case that the target is located at a hub node, a walker 
performing MERW arrives at the destination more quickly than 
performing TURW. However, for the two cases that target node is 
a node with the smallest degree or a node selected uniformly, MERW 
is less efficient for finding the target than for TURW. 

We have also found that the values of MFPT for MERW in scale- 
free networks are distributed over a larger range than their counter- 
parts for TURW. Thus, as an indicator of node importance, MFPT 
for MERW is better at discriminating influential nodes from com- 
mon noncentral nodes. Finally, we note that our approximate ana- 
lytical results only hold for uncorrelated networks. Since in real 
networks, there exist ubiquitous degree correlations among nodes^", 
it would be interesting to extend our methods to correlated networks 
in the future. 

Methods 

Expressing FPT for MERW in a network in terms of the spectra of its adjacency 
matrix. It has been reported^*" that for generic random walks in a weighted network, 
the FPT Tjj from node i to node j can be represented by the eigenvalues and 
eigenvectors of the following matrix F defmed as 



r=s 5WS2=S5PS 5, 



(32) 



where S is the diagonal strength matrix with its ith diagonal entry equal to the strength 
s, of node i. It is evident that matrix F is real and similar to the transition matrix P and 
thus has the same set of eigenvalues as P. Let ^^,^2 , ■ ■ ■ ,^3 be the N eigenvalues of 
matrix F for a network of size N, rearranged as 1 — > > ■ ■ ■ >^lj? and let 
i/'i,t//2, ■ ■ ■ .yjy^ denote the corresponding normalized and mutually orthogonal 
eigenvectors, where ij/^ — (i/'ji,t//,25 ■ ■ ■ ^tAiw)^- Then, the FPT for a walker starting 
from node / to first arrive at node j can be expressed as^^ 



(33) 



For the particular weighted network associated with MERW, its weight matrix W 
satisfies 



W= — SiASi. 



Inserting equation (34) into equation (32) leads to 



(34) 



(35) 



which implies that the eigenvalues and eigenvectors of the two matrices F and A, 
satisfy the following one-to-one relations: 



and 



\j/i = Hi. 



(36) 



(37) 



kiP{kj\ki) 



(41) 



where P{kj\ki) is the conditional probability^^ that a node of degree kj is directly 
connected to a node with degree kj. Note that A can be also interpreted as a weight 
matrix of a weighted, fully connected graph, which is obtained by annealed network 
approach^". Since fi is the eigenvector of A corresponding to the principle eigenvalue 
Ai, we can approximate /ti by J^i. Next, we evaluate the principle eigenvector Ji of A. 

Since in an uncorrelated network, the degrees of the two nodes connecting any edge 
are completely independent, the conditional probability can be estimated as 



p{kj\k>)=kjp{kj)/{d). 



(42) 



where {d) is the average node degree. Instituting equation (42) into equation (41) 
yields 



"kikj - 



kikj 



which indicates that the matrix A can be represented as 



A=— -kk', 

N{d) 



(43) 



(44) 



where k is the degree sequence of the network and can be denoted by a vector as 

k = (fci,*:2,...,/:„)\ (45) 

For a matrix having the form txa ^ , where a is a nonzero vector, its rank is 1 . Therefore, 
the rank of A is 1, and A has exactly one nonzero eigenvalue 



A = tr(A) = — I— V/t?. 



(46) 



Having obtained the principal eigenvalue we continue to determine its corres- 
ponding eigenvector //. Obviously, A and Ji satisfy the following relation: 



which can be reexpressed as a system of equations: 

' kiji^+kijii^ ^knji^ = -^^Jii, 

kipi+k2p2 H hfcw/l„ = ^Ji2, 

^ kifii+k2fi2+ h = j^Aiw, 

where ^i, is the ith entry of /i. Therefore, 



Ih 



kN ' 



(47) 



(48) 



(49) 



Combining with normalized condition — 1> equation (49) can be solved to 

obtain 



(50) 



Approximate by ^i^ leads to 



Substituting equations (36) and (37) into equation (33) and considering 5/ — Ai/f,;, we 
obtain 



^1} sr-^N 



(51) 



_ 1 ^ ;.i 



^'kj ~ l^kil-'kj 



(38) 



which provides a close-form expression of the FPT for MERW starting from an 
arbitrary node /' to another node j. 

Approximation of the principal eigenvector. Mean-field theory assumes that nodes 
having the same degree share the same structural properties'^. Based on this 
hypothesis, we use /((fc) to denote the value of elements of principal eigenvector 
corresponding to nodes with degree k. By definition. 



A/(i-Ai^(i. 



(39) 



Applying the coarse-graining idea to degree classes'^, equation (39) is equivalent to 

Am-;.i/1, (40) 

where fX — (i.i(ki),fi(k2), ■ ■ ■ ^ii{ki^)Y . The entry 0.^.^. of matrix A defines the 
probability that two nodes of degree /c, and k^ are adjacent, namely 
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